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ABSTRACT 

Wide binaries made up of two white dwarfs (WDs) receive far less attention than their tight coun- 
terparts. However, our tests using the binary population synthesis code Star Track indicate that, for 
any set of reasonable initial conditions, there exists a significant observable population of double white 
dwarfs (WD WDs) with orbital separations of 10 2 to 10 5 AU. We adapt the technique of Dhital et al. to 
search for candidate common proper motion WD companions separated by <10' around the >12,000 
spectroscopically confirmed hydrogen-atmosphere WDs recently identified in the Sloan Digital Sky 
Survey. Using two techniques to separate random alignments from high-confidence pairs, we find nine 
new high-probability wide WDWDs and confirm three previously identified candidate wide WDWDs. 
This brings the number of known wide WDWDs to 45; our new pairs are a significant addition to 
the sample, especially at small proper motions (<200 mas/yr) and large angular separations (>10"). 
Spectroscopic follow-up and an extension of this method to a larger, photometrically selected set of 
SDSS WDs may eventually produce a large enough dataset for WDWDs to realize their full potential 
as testbeds for theories of stellar evolution. 
Subject headings: binaries: general — white dwarfs 



1. INTRODUCTION 

Hydrogen-atmosphere (DA) white dwarfs (WDs), the 
evolutionary endpoints for stars with main-sequence 
masses between 0.8 and 8 Mq, are second in num- 
ber only to lo w-mass main-sequence s tars in the Solar 
neighborhood ( Rowell fc Hamblvll2011l ). However, while 
they can remain at relatively high temperatures for Gyr, 
WDs are small and therefore usually faint objects, so 
that assembling complete catalogs of WDs has histor- 
ically been challenging. The last decade has seen an 
impressive growth in our ability to find and character- 
ize WDs, l argely thanks to the Sloan Digital Sky Sur- 
vey (SDSS; lYork et al.l [2 000l). The first SPSS WD cata- 
log included 2551 WDs (|Kleinman et al.ll2004h ; the sec- 
ond nearly quadrupled that number, a nd included 6000 
new, spectroscopically confirmed WDs (Eisenstei n et al.l 
12006ft . The most recent SPS S WD catalog, based on the 
SDSS Pata Release 7 fPR7: IAbazaiian et al.ll2009f ). con- 
tains ^20,000 stars, of which 12,538 arc classified as PAs 
(Kleinman et al., in prep.). 

Unsurprisingly, many of these WPs are in binary 
systems. The zoo of SPSS WP binaries is diverse; 
among its largest pens are those containing spectro- 
scopic WP/M-dwarf pairs (e . g.. ISilvestri et al.l [20071 ; 
iRebassa-Mansergas et alj[2010t Morgan et al.l 20121) and 
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ISzkodv et al.l 1201 If ). A smaller but especially inter- 
esting pen is the one for double WP (WPWP) sys- 
tems. Pepending on the composition and mass of 
the two WPs, WPWPs may be a formation chan- 
nel for AM CVn-type bin aries (also found in SPSS; 
I Anderson et alJ 120051 . [2008) , millisecond pulsars (with 
the neutron stars for med through a ccretion-induced col- 
lapse; cf. discussion in lNelson1l2005D . and Type la super- 
novae (in the double-d e generate scenario; lWebbinldll~98l 
llben fc Tutukovl Il98l Ivan Kerkwiik et alJ |2010T T Fur- 
thermore, the discovery of WPWPs with mass ratios far 
from unity and with at least one WP with M < 0.2 Mq 
tests our understanding of common envelo pe evolution 
(jAgiieros et al.ll200l iKilic et alj|2010l l201lh . 

Most known SPSS WPWPs are unresolved and/or 
compact binaries. By comparison, the expected popu- 
lation of wide, resolved WPWPs remains mostly unex- 
plored, largely because of the observational challenges 
these types of binaries present, whatever their compo- 
nent stars. Wide binaries can have orbital periods that 
are much longer than human lifetimes. Members of these 
systems can only be identified astrometrically and con- 
firmed through radial velocity measurements, although 
in practice the latter are hard to obtain for large samples 
of field stars Q Over his 70 year career, Luyten perfected 
a method that relies on proper motion (//) measurements 
to identify binar y systems (see for example the discussion 
in lLuvtenlll988f ). The resulting pairs are known as com- 
mon proper- motion binaries (CPMBs). This method is 
now commonly used and generally relies on high-proper- 
motion, magnitude-limited catalogs, s o that searches for 
CPMBs tend to uncover nearby stars (jPhital et al.ll2010L 
and references therein). 

7 Furthermore, typical line-of-sight velocities for WDs are com- 
parable to the gravitational re dshift experienced by the photons 
they emit jSilvestri et al.ll200U ). 
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Wide binaries provide a unique perspective on stellar 
evolution: for those with orbital separations a > 10 2 
AU, the coeval components are far enough apart that 
mass exchanges are unlikely to h ave significantly im- 
pacted their indiv idual evolutions (|Silvestri et alj [20011 
IFarihi et al.l f2006). Wide binaries that include a main- 
sequence star and a WD can be used to determine 
the initial-to-final mass relation (IFM R) for WDs (e.g., 
iCatalan et all 120081 : iZhao et al.ll2012fl . The WD mass 
and cooling age are determined spectroscopically, and 
this age is subtracted from the main-sequence star's age 
to determine the WD's main-sequence lifetime and hence 
initial mass. WD WDs, for which total ages are not gener- 
ally availabl e, can still provide robus t constraints on the 
IFMR (e.g.. iFinlev fe Koesteil Il99l . WDWD systems 
can also be used to constrain the effects of mass-loss on 
orbital evolution, for example by comparing their sep- 
arations to the distributions for binaries containing two 
main -sequence stars a nd a main-sequence star and a WD 
(e.g.. lSion et al.|[l991l ). Unfortunately, the current sam- 
ple of wide WDWD systems is small and heterogeneous, 
severely limiting i ts utility for t hese s tudies. 

Th e surveys of iGiclas etail (j!971[ ) and iLuvtenl (e.g., 
1979) identified several thousand CPMBs; of these, a 
few hundred were thou ght to contain a WD. Spectro- 
scopic follow-up allowed | Greensteinl ([1986) to confirm six 
as WD WDs and lSion et alj (|1991f ) to increase that num- 
ber to 21. Since then, individual WD WDs have been 
rep orted (e.g., those found in the Palomar Green sur- 
vey; ^rihreOI][200^) , but there has been no systematic 
search for new wide WD WDs. 

SDSS is an excellent dataset for identifying wide WD- 
WDs. The photometric survey covered a large area 
to unprecedented de pth (>10 deg 2 and ~22 mag; 
lAbazaiian et al.ll2009D , and matches between SDSS and 
USN O-B have been use d to generate a proper motion cat- 
alog (|Munn et al.ll2004D that is i ntegr ated into the SDSS 
database. Indeed. iDhital et alj (|2010D used SDSS to un- 
cover wide companions to low-mass main-sequence stars 
based on common proper motions: these authors identi- 
fied over 1300 CPMBs, of which 21 i nclude one WD. 

We adapt the method developed bv lDhital et al.l (|2010l ) 
to search for widely separated WD companions to the 
set of spectroscopically confirmed DA WDs included in 
the forthcoming Kleinman et al. catalog. In Section 2 
we discuss the population synthesis predictions for the 
orbital distribution of WD binaries that motivated this 
search. We present in Section 3 our method for identi- 
fying candidate CPMBs and in Section 4 the properties 
of our newly discovered WD WDs; we also compare these 
properties to those of the previously known pairs. We 
conclude in Section 5. 

2. MOTIVATION 

Whether a WDWD is observed today as a tight or 
a wide binary depends primarily on the orbital separa- 
tion of the binary at birth. Because WD progenitors go 
through mass-losing giant phases, WDWD orbits may 
expand up to a f actor of five due t o conservation of angu- 
lar momentum (fGrecnstcin 1986). If, however, the WD 
progenitors overflow their Roche lobes while in a giant 
phase, they will likely enter unstable mass transfer, caus- 
ing a rapid, order-of- ma gnitude shrinking of the orbit 
(e.g.. Ilben fe Liviol 119931 ). We therefore expect to ob- 



TABLE 1 

Initial Conditions for Population Synthesis 



Parameter 



Distribution 



Range 



Eccentricity 

Primary IMF 

Mass ratio 
Orbital separation 



#i oc 2e 

ae 

dM * \ 
dq 

IT oc a- 1 

da 



M 
M- 



-2.2 



< e < 1 
0-5 < *5 < 1 

< q < 1 
R L <a< 10 5 AU 



Note. — Rl is the radius of the Roche lobe. 

serve a bimodal distribution of orbital separations, with 
a population of wide WDWDs that avoided any mass 
transfer phases and one of tight WDWDs that underwent 
unstable mass transfer. We use population synthesis to 
characterize this distribution in detail. 

2.1. Population Synthesis of WDWD Systems 

Population synthesis is commonly used to analyze 
the formation and evolut ion of WDWDs. For exam - 
ple. iNelemans et al.l (|200lD and lNelemans fc Tout! (|2005T ) 
have used it to ex amin e the physic s of m ass transfer, 
IRuiter et all (|2011D and IMeng et al.l (|2011Q to estimate 
SN la rate s due to WDWD m ergers, and IRuiter et al.1 
(j2010l ) and lYu fc Jeffervl (|2010D to predict the gravita- 
tional wave signal from inspiraling WDWDs — all phe- 
nomena related to tight WDWDs. Wide WDWDs, how- 
ever, have largely been ignored in these studies. 

Using the binary populatio n synthesis code StarTrack 
(|Belczvnski et al.1 l2002. 2008j), we evolve a sample of 10 6 
zero age main sequence (ZAMS) binaries and determine 
the orbital separation distribution of the ^5% that be- 
come WDWDs. We summarize our initial conditions in 
Tabic Q] and describe them briefly below. 
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Fig. 1. — Normalized orbital separation distribution of binaries 
produced by our population synthesis code once they become WD- 
WDs (solid line) and after perturbing effects have been taken into 
account (dashed line). The compact systems are depleted due to 
mergers, while the widest systems are disrupted due to three-body 
interactions. 

We use t he standard Ambarts u mian eccentricity dis- 
tribu tion (|Ambartsumianl 119371 : iDuauennov fc Mavorl 
I1991D . The initial mass of the more massive primary 
star is drawn from a Kroupa initial mass function (IMF; 
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iKroupa et alj fl993T) . and the initial mass of the com- 
panion star is a random fraction of the primary mass 
(|Kobulnickv fe Frverj|2007D . Because of the weak depen- 
dence of the orbital separation on the component mass 
and orbital eccentricity, we expect that other (reason- 
able) distributions for these variables would not signifi- 
cantly alter the orbital separation distribution. 

The initial orbital separations range from the two 
stars starting just outside of contact to being 10 5 AU 
apart. The maximum separation of binaries at birth is 
not well constrained. However, ~15% of all G dwarfs 
are found in binary system s with separations >10 AU 
(|Duauennov fc Mavorll99l1), the mean radii of pre-stellar 
cores are ~10 5 AU ([Clemens et al.l fl99l . and other 
studies have ident ified binaries with a ~ 10 5 AU (e.g., 
IDhital et all2010D , so that this maximum initial separa- 
tion is a reasonable estimate. 

We choose a logarithmical l y flat distribution for 
a (|Qpikl 119241 : iPoveda et ail I2007T) . Observation- 
ally, the distribution of bi rth orbital separations may 



be very different (e.g. 
IChaname fc Gouldl l20f)1 



Duauennov fc Mayor! [1991; 
Lepine fc Bongiorna 2007). 
Howev er, when we test the broke n power-law distribu- 
tion of iLepine fc Bongiornol (|2007h . we see no significant 
differences in the resulting present day orbital separation 
distribution of WDWDs. 

The solid line in Figure [1] is the orbital separation dis- 
tribution immediately after the birth of the second WD. 
As expected, this distribution is bimodal, with the num- 
ber of systems with a > 10 2 AU dominating the over- 
all population. Interestingly, about 10% of these wide 
pairs did experience mass transfer: the primary under- 
went stable Roche lobe overflow while on the asymptotic 
giant branch. In 90% of cases, by contrast, the two stars 
can be considered to have evolved independently. 

While the number and distribution of tight systems in 
Figure [1] is strongly dependent upon e.g., the prescrip- 
tion used to describe common envelope evolution, the 
existence of a population of a population of WDWDs 
widely separated at birth is a robust prediction of our 
code. 

2.2. Evolution of the Synthesized WDWD Orbits 

WD WD orbits are modified over time by gravitational 
wave radiation and weak interactions with other bod ies 
in the Galaxy. We use the equations of lPeters! (|1964l ) to 
model the first of these effects, which is most important 
for the tightest binaries in Figure [T] We model the sec- 
ond effect using the Fokker-Planck approximation, which 
determines the diffusion of energy into a nd out of the bi- 
nary s ystem. Using this approximation, Weinber g et al.l 
( 1987) found the characteristic lifetime t of a binary con- 
sisting of two 0.6 Mq WDs to be: 



t 



2.9(—^—) 'Gyr. 
V10 5 AU/ y 



(1) 



The lifetime of a canonical WDWD with a ~ 10 5 AU 
is therefore a few GyrQ 

8 This may be an overestimate of a wide binary's lifetime be- 
cause it ignores the dif ferential pull of the Galactic potential 
llJiang Trcmainc 2010). Although this effect is important, ap- 
plying these authors' semi-analytic method to our population syn- 
thesis output is beyond the scope of this paper. 



We assign each ZAMS binary a birth time uniformly 
distributed over the lifetime of the Galactic disk (0-10 
Gyr). The dashed line in Figure [T] shows the distribution 
of surviving binaries at the present day (10 Gyr), af- 
ter gravitational radiation and Galactic interactions have 
caused some to merge and disrupted others. Even with 
our gross overestimate of the effect of Galactic perturba- 
tions, a large number of the wide WDWDs survive and 
should be observable today. 
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Fig. 2. — Expected cumulative percentage of WDWD as a func- 
tion of separation once projection effects are taken into account. 
The binaries are artificially placed at a distance of 250 pc. Close 
pairs dominate the distribution, but ~15% have 8 > 2'. 



2.3. Predicted Observed Angular Separation Distribution 

To translate our synthesized population into an ob- 
servable population, we assign each WDWD an argu- 
ment of pericenter, mean anomaly, and cosine of the in- 
clination angle, all randomly selected from flat distribu- 
tions. These parameters, combined with the intrinsic 
orbital parameters of the binary (eccentricity, masses, 
and orbital separation) define the orbital positions of the 
two WDs. Using a fiducial distance of 250 pc, we then 
determine the projected angular separation (8) of each 
WDWD. Binaries with 8 < 8" or >10' are eliminated 
because our method for identifying WDWDs (described 
below) is not sensitive to these. The predicted observed 8 
distribution of our synthesized WDWDs is shown in Fig- 
ure [21 We find that the qualitative characteristics of this 
distribution are conserved for any reasonable assumption 
for the distance to the binaries. Although the WDWD 
population is dominated by pairs with separations <1', a 
small fraction is expected to be observed at separations 
extending to 10'. 

3. IDENTIFYING WIDE WDWDS IN SDSS 

3.1. Candidate Binary Selection Process 

The SDSS data reduction pipeline matches objects to 
objects in the USNO-B catalog, which has a limiting 
magnitude of ~21 mag. Proper motions calculated dur- 
ing this matching are included in the propermotion table 
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Fig. 3. — (g — r) vs. (u — g) for 8605 spectroscopically confirmed 
DAs with a unique entry in the new Kleinman et al. catalog and 
with ugriz photometric errors <0.15 mag (stars and dark con- 
tours). The points and light-colored contours to the upper right 
indicate the colors of >99,000 randomly selected stars with ugriz 
errors <0.1 mag, ugriz > 15.5 mag, and g < 20 mag (these stars 
also meet the same proper motion constraints that were applied 
to our primary WDs). These stars are included to highlight the 
small overlap between the WD locus and the main sequence. Over- 
plotted are different regions used to photometrically s elect WDs 
and A stars. The region denned by Girvcn ct al. ( 2011) traces the 
empirical DA locus extremely well and includes 96% of the 8605 
DAs, a significantly l arger fraction than returned by the standard 
IRichards et al. (2002) color-cuts. 

of the DR7 catalog^ We use the CasJobs databasj^l to 
obtain all the available proper motion data for the DA 
WDs in the Kleinman et al. catalog: 11,563 DAs (92%) 
have measu red proper m o tions. 

Following iMunn et alj (|2004D , we eliminate any WD 
with a rms fit residual <350 mas in either right ascen- 
sion (a) or declination (5) and any WD with more than 
one possible USNO-B counterpart (i.e., match ^ 1). We 
also eliminate any WD with a total /i < 20 mas/yr or 
with > 10 mas/yr in either coordinate. ~56% of 
the DAs in the Kleinman et al. catalog do not survive 
these cuts; most of these WDs have total proper motions 
below our threshold. We refer to the remaining ~5500 
DAs as primaries in our CPMBs and to their candidate 
companions as secondaries. 

Next, we query CasJobs for all objects offset by <10' 
from each of our primary WDs, yielding ^4.1 x 10 6 ob- 
jectsE3 Applying the proper motion quality cuts de- 
scribed above pares our list of potential secondaries down 
to ^2.2 x 10 5 objects. (Most of the eliminated objects 
lack proper motions only because they are too faint to be 
included in the USNO-B catalog.) We eliminate objects 
with poor photometry, i.e., with g-band errors >0.1 mag 
or uriz errors >1.0 mag; this leaves ~2.0 x 10 5 objects. 

We then search for objects with a proper motion 
matching that of our primary WDs. Typical proper mo- 
tion observational errors are ~4 mas/yr. In the most 

9 Proper motions generated from this matching before DR7 
contained a syst ematic error that has since been corrected 
HMunn et al.ll200"8f ) . 

10 http://cas.sdss.org/casjobs/ 

11 We require that these objects be defined by SDSS as primary, 
m eaning they met a num ber of quality criteria that are described 
in lStoughton et al.l ((2002T ). 



TABLE 2 

Candidate Wide WDWDs with SDSS Spectra for Both 
Components 







Di 


D 2 


RVi 


RV 2 


Name 


(pc) 


(PC) 


(km/s) 


(km/s) 


J0332- 


-0049 


143 ± 21 


194 ± 29 


-31 ± 11 


-45 ±7 


J0915- 


-0947 


216 ± 32 


115 ± 17 


-58 ±8 


±18 ± 7 


J101H 


-2450 


436 ± 65 


772 ±116 


±54 ±17 


±9 ±30 


J1113H 


-3238 


147 ± 22 


106 ± 16 


±63 ± 36 


±26 ± 35 


J1257n 


-1925 


516 ± 77 


457 ± 69 


±23 ± 27 


±33 ± 22 


J1309n 


-5503 


137 ±21 


83 ±12 


-50 ± 23 


±45 ± 15 


J1555- 


-0239 


291 ± 44 


143 ± 21 


±11 ± 14 


-3 ±27 


J2326- 


-0023 


115 ± 17 


101 ± 15 


-23 ±62 


-11 ± 9 



Note. — These are all confir med DA WDs included in 
either the [Eiscnstein ct al. (2006) or Kleinman et al. SDSS 
catalogs. 



extreme case of a face-on orbit at the observational limit 
9 = 8" and at a distance of 50 pc, two WDs with typical 
M = 0.6 Mq in a circular binary will have an orbital 
velocity of ~0.8 km/s. This translates to a differential 
fi ~ 3.4 mas/yr, which is of the same order as our obser- 
vational (7 U - Since all our binaries are found at distances 



>50 pc (see below) and have 9 > 
of our binaries can be ignored. 



the orbital velocities 



We define a proper motion match in the same manner 



we aenne a proper n 
as IDhital et all (2010): 



&Au a 



( A(i s 

\. a Au s 



< 2, 



(2) 



where A^i is the scalar proper motion difference in a 
and £, and is the error in the corresponding A/i, 
derived from the quadrature sum of the individual er- 
rors in fi. This greatly reduces the number of candidate 
secondaries, to only 7129. 

Finally, to identify candidate WDs among these sec- 
ondaries, we select objects within 0.5er of the region 
in (q — r) yersus ( u — g) color-color space defined by 
iGirven et all (|2011h as occupied by WDs. This region 
encompa sses a larger area of t he WD locus than do the 
standard Rich ards et alJ (|2002l ) color cuts. Furthermore, 
the lGirven et alj (|2011l ) region does not o verlap with the 
(a — r ) ver sus (u — g) b oxes defined by IRichards et al.l 
(|20"02T) and lCovev etaTI ((2001 for main-sequence A and 
early F stars, the most likely stellar contaminants (see 
Figure EDB 41 SDSS objects survive this final cut. 



3.2. Distances (and Radial Velocities) 

In addition to matching proper motions, bound bina- 
ries should have matching distances and (when avail- 
able) radial velocities (RVs). For our candidate pri- 
maries, we use the spectroscopically derived T c ff and 
log g values from the Kleinman et al. catalog and lin- 
early interpolate T e g and quadra tically interpolate log 
q in model evolutionary grida^l (T remblav et alJ |2011l : 
iHolberg fc Bergeronll2006HKowalski fc Saumonll2006l ) to 
determine their absolute magnitudes in each of the five 
SDSS bands, taking into account the best fit for Galactic 

12 The separation between these regions is cleanest in this com- 
bination of colors. 

1 ' i http: //www. astro .umontreal . ca/~bergeron/CoolingModels/ 
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Fig. 4. — 8 distribution of real (solid line) and of false (dashed) 
candidate CPMBs. The primary WD positions were shifted by 
±1° in both a and <5; we then applied our method for identifying 
common proper motion companions to these shifted stars. There 
is an excess in the distribution of real candidates for 8 < 1.5 — 2'; 
at larger 8, the contamination by false pairs is essentially 100%. 

rcddening'0 The distance to each primary is then the av- 
erage of the distances derived in each band [3 Formally, 
the uncertainties associated with these distances are very 
small (a few percent). However, these do not reflect the 
hard-to-quantify uncertainties in the model grids, and we 
therefore adopt a more realistic uncertainty of 15% for 
these distances. 

Determining photometric distances to our candidate 
secondaries, which generally lack SDSS spectroscopy, is 
less straightforward. While T e g can be derived from 
SDSS photometry to within a few 100 K, the typical 
photometric errors are sufficient to create uncertainties 
of order half a decade in the derived log g. As a result, 
one typically assumes log g — 8.0 (the Kleinman et al. 
DA sample is strongly peaked at this value). We compare 
the distances calculated in this manner for our candidate 
primaries to those derived from their spectroscopically 
determined T c ff and log g values. We find that they are 
consistent only to within a factor of 2. In general, there- 
fore, the photometric distance uncertainties are too large 
for us to differentiate between real and spurious binaries 
based on comparisons of the distance to the two stars. 

There are eight candidate pairs in which both compo- 
nents have SDSS spectra. We provide the spectroscopi- 
cally derived distances and cataloged RV measurements 
for these WDs in Table [2] and discuss these systems in 
more detail in Section 0J 

3.3. Estimating the Purity of Our Sample 

We use two methods to test the robustness of our selec- 
tion. First, we estimate the overall contamination of our 
sample by shifting the positions of our primary WDs and 
using the same photometric and proper motion criteria 
to identify (false) companions to these shifted primaries. 
Second, we use the Monte Carlo approach developed by 

14 One candidate primary's T a g and log g values in the Kleinman 
ct al. catalog were unusable because the code we use to generate 
these values rejected the associated calibrations. We obtain T c ff 
and log g estimates for this WD by redoing the fit (K. Oliveira, 
private communication). 

15 Fits to the grids also provide cooling ages and masses. 
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Fig. 5. — 8 distribution of the 41 actual WDWD candidates 
(solid line) compared to that of the artificially generated candidate 
CPMBs (dashed line). The dotted line is the predicted distribution 
from population synthesis (shown in Figure [2} normalized to the 
first bin in the real distribution, and suggests that our sample is 
incomplete at 8 < 2'. Our distribution of real candidates includes 
more candidates at large 8 than is expected from our population 
synthesis results, and is strongly contaminated for 8 > 2'. 

iDhital et al.l ((2010) to estimate how likely it is to find 
a random star whose characteristics match those of our 
candidate primaries. 

3.3.1. Empirical False Positive Determination 

We shift the positions of our primary DA WDs four 
times (±1° in both a and 5) and use the [i cuts described 
above to identify candidate proper-motion companions 
to these stars. In Figure 2J we compare the normalized 
9 distribution of the resulting population of false can- 
didate CPMBs after these shifts and that of the 7129 
candidate CPMBs identified when using the true DA po- 
sitions. There is an excess in the distribution of real 
candidates for 9 < 1.5 — 2', while at larger separations, 
the distributions are equivalent, suggesting that most of 
our actual candidates are due to random matches. 

In Figure we compare the 9 distribution for the can- 
didate WDWDs selected from the real CPMB candidates 
to that for the false candidate CPMBs shown in FigurelH 
Given the small number of candidate WDWDs, we also 
plot the predicted distribution for WDWDs from Fig- 
ure [5] normalized to the first 9 bin in our distribution of 
real candidates. While Figure 0] suggests that true bina- 
ries are most likely at separations <1', the total expected 
number of observed wide WDWDs (dotted line in Fig- 
ure [5]) is significantly larger than the expected number of 
false positives (dashed line) out to 9 ~ 2'. Although none 
of our actual candidates have 1' < 9 < 2', our results sug- 
gest that such WDWDs, if found in future searches using 
similar constraints, are likely to be true binaries. 

3.3.2. Galactic Model 

Each realization of the IDhital et al.l (|2010l ) Galactic 
model populates a 30' x 30' conical volume centered at 
the position of the primary WD up to distances of 2500 
pc from the Sun. The model assigns a position in six- 
dimensional phase space to each simulated star, assuming 
three kinematic components of the Galaxy corresponding 
to thin disk, thick disk, and halo populations. 
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Fig. 6. — Percentage of model realizations that return a random 
star with the same properties as those of the 41 primary WDs in 
our candidate WDWD binaries; 6 is the separation of the corre- 
sponding candidate WDWDs. The arrows indicate WDs for which 
P m is either less than or greater than the y-axis limits. P m can 
be >100% because a given realization may have more than one 
star with properties matching those of the candidate being tested. 
Note that there is a primary at 6 ~ 5' with P m ~ 100% and two 
primaries at ~ 8' that are also partially covered by the arrows. 
13 primaries have P m < 1% and 9 < 2'. 

After 10 5 realizations of the model, we count the num- 
ber of rendered stars for which /i matches fi pr im,ary (as 
defined by Equation^ and that are at a separation from 
the primary smaller than or equal to the separation of the 
corresponding candidate binary. A matching star also 
has to be at a distance consistent with the distance to the 
primary WD to within the quadratic sum of the distance 
uncertainties to the primary and secondary in the binary 
being tested. In most cases, this corresponds to search- 
ing a relatively large volume along the linc-of-sight for 
matches, since while the spectroscopic distances uncer- 
tainties are 15%, the photometric distance uncertainties 
are taken to be 100%. For the eight primaries whose sec- 
ondaries also have spectroscopic distances and RV mea- 
surements (see Tabled, the searched volume is smaller, 
and we further require that a match have an RV consis- 
tent with that of the primary WD to within la. Unsur- 
prisingly, this results in a systematically lower number of 
random matches to the primaries in these pairs. 

Figure [B] shows the percentage P m of realizations that 
return a random star whose properties match those of our 
41 candidate primary WDs. These results are in agree- 
ment with those from our first test: finding a random 
match within 2' of one of our primary WDs is extremely 
unlikely, while pairs with 9 > 5' are much more likely 
to be random matches. Interestingly, however, this test 
suggests that there are systems with 2' < 9 < 5' that are 
more likely to be real than random, as the likelihood of 
finding a random star with properties matching those of 
the WD primaries is small. 

4. RESULTS AND DISCUSSION 
4.1. New WDWDs 

Based on the tests described above, we set 9 < 2' and 
P m < 1% as our criteria for identifying true pairs. 

Of the candidate pairs with SDSS spectra for both 
WDs, J101 1+2450 and J1555+0239 have 9 > 2' and 




-10 12 3 

(g4) 

Fig. 7. — Reduced proper motion as a function of (g — i) for the 
SDSS stars presented in Figure \3\ (points and contours) and for 
members of our high-confidence candidate WDWDs. Spcctroscop- 
ically confirmed WDs are shown as blue stars, while the red stars 
lack spectra. The Vt = 30 km/s line marks the expected loca- 
tion of disk WDs and the Vt = 150 km/s line that of halo WDs 
(Kilic ct al. 2006). All of our candidates lie in the expected region 
for WDs, far from the main-sequence locus. 

Pm > 1%, and are therefore very likely to be random 
alignments. The distance and RV data in Table [2] are 
consistent with this, as the disagreements in one or both 
measurements are significant. J0915+0947, for which 
9 = 4.5', does have P m < 1%, but the RV measure- 
ments for the two WDs are highly discrepant (—58 ± 8 
and +18±7 km/s), implying that this is a random match. 

Conversely, five pairs listed in Table [2] have 9 < 2' 
and P m < 1%: J0332-0049, J1113+3238, J1257+1925, 
J1309+5503, and J2326-0023. We check the distances to 
the two components and the RV measurements for these 
for consistency: of the five, J1309+5503 is the only one 
for which the disagreements between the measurements 
for each DA (— 50±23 and +45±15 km/s) are significant 
enough to eliminate it from our list of candidate pairs. 
We identify the other four pairs as high-confidence can- 
didate WDWD systems and present them along with the 
eight systems that lack SDSS spectra for the secondaries 
but also meet our criteria in Table [3] 

In addition to the four pairs with SDSS spectra for 
both DAs, a SIMBAE0 search finds that two sys- 
te ms in Table [3] have s econdaries classified as DAs 
by iMcCook fc Sionl (|1999f ) . Spectroscopic follow-up is 
clearly needed to confirm the nature of the secondary 
in th e remaining s i x syst ems. For objects 17 < g < 19 
mag, iGirven et al.l (|2011h estimate that their photomet- 
ric selection leads to a quasar contamination rate of 
<35%, primarily at the faint end. At the bright end, 
IGirven et alj (|201lD find that early-type main-sequence 
stars and subdwarfs are the main contaminants, and es- 
timate their contamination rate to be <20%. 

However, the contamination rate drops sharply when 
a proper motion constraint is applied: if one imposes 
[i > 20 mas/yr, the quasar contamination rate becomes 
negligible and ~90% of the hot stars are eliminated. Con- 
sidering that our high-confidence pairs all have fj, > 30 
mas/yr, we expect that no more than one of our high- 

16 http://simbad.u-strasbg.fr/simbad/ 
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TABLE 3 

High-Confidence Candidate Wide WDWDs 



Name 


a 


6 


9 (mag) 


fi a (mas/yr) 


fj, s (mas/yr) 


Type 


6 (") 


1 m 


Ref. 


J0000-1051 


00:00:22.5 
00:00:22.8 


-10:51:42.1 
-10:51:26.6 


18.91 ± 0.04 
20.21 ± 0.04 


+45.3 + 4.6 
+42.3 ± 4.1 


-25.3 + 4.6 
-24.3 + 4.1 


DA 


16.1 


<0.1% 


1,2 


J0029+0015 a 


00:29:25.6 
00:29:25.3 


+00:15:52.7 
+00:15:59.8 


18.48 ± 0.03 
19.59 ± 0.02 


-28.9 + 3.1 
-27.9 + 3.6 


-23.0 + 3.1 
-23.9 ± 3.6 


DA 


8.6 


<0.1% 


1,2 


J0332— 0049 


03:32:36.6 
03:32:36.9 


-00:49:18.4 
-00:49:36.9 


18.20 + 0.03 
15.64 + 0.02 


-24.8 + 5.4 
-30.9 + 4.5 


-23.6 + 5.4 
-23.3 + 4.5 


DA 
DA 


18.9 


<0.1% 


1,2,3,4 
2 


J1002+3606 


10:02:44.9 
10:02:45.8 


+36:06:29.5 
+36:06:53.3 


18.92 ± 0.03 
19.04 + 0.02 


-32.9 + 3.3 
-29.9 + 3.4 


-27.5 ± 3.3 
-27.0 + 3.4 


DA 


26.5 


0.12% 


1,2 


I1054+5307 b 


10:54:49.9 
10:54:49.2 


+53:07:59.2 
+53:07:15.2 


17.92 ± 0.04 
17.52 ± 0.03 


-113.9 + 3.1 
-112.9 ± 2.9 


-38.4 + 3.1 
-36.0 + 2.9 


DA 
DA 


44.5 


0.12% 


1,2,4,5 
4,5 


J1113+3238 c 


11:13:19.4 
11:13:22.6 


+32:38:17.9 
+32:38:58.9 


19.03 ± 0.03 
19.12 ± 0.04 


-162.9 + 3.1 
-158.8 ± 3.2 


+58.0 + 3.1 
+58.0 + 3.2 


DA 
DA 


56.7 


<0.1% 


2 
2 


J1203+4948 


12:03:11.5 
12:03:11.0 


+49:48:32.4 
+49:48:50.8 


19.03 ± 0.03 
17.35 ± 0.02 


-97.6 + 3.4 
-98.2 + 2.9 


-36.5 ± 3.4 
-39.1 ± 2.9 


DA 


19.3 


<0.1% 


2 


J1257+1925 


12:57:20.9 
12:57:21.1 


+19:25:03.7 
+19:24:51.8 


19.88 ± 0.06 
17.07 ± 0.03 


-38.4 + 5.4 
-33.0 ± 2.7 


-31.6 + 5.4 
-33.5 ± 2.7 


DA 
DA 


12.2 


<0.1% 


2 
2 


J1412+4216 d 


14:12:08.9 
14:12:07.7 


+42:16:24.6 
+42:16:27.1 


18.46 ± 0.02 
15.83 + 0.01 


-80.3 + 3.1 
-81.7 + 2.7 


-57.5 ± 3.1 
-61.4 + 2.7 


DA 
DA 


13.6 


<0.1% 


1,2 
4,6 


J1703+3304 


17:03:55.9 
17:03:56.9 


+33:04:38.3 
+33:04:35.8 


18.81 ± 0.02 
18.16 + 0.01 


-1.8 + 3.4 
+0.3 + 3.1 


-51.2 ± 3.4 
-50.5 + 3.1 


DA 


11.2 


<0.1% 


1,2 


J2115-0741 c 


21:15:07.4 
21:15:07.4 


-07:41:51.5 
-07:41:34.5 


17.47 ± 0.02 
16.81 ± 0.01 


-25.2 + 2.9 
-30.0 + 2.8 


-117.2 + 2.9 
-117.9 + 2.8 


DA 


17.0 


<0.1% 


1,2 


J2326-0023 6 


23:26:58.8 
23:26:59.3 


-00:23:39.9 
-00:23:48.1 


19.33 ± 0.05 
17.49 ± 0.02 


+51.6 + 3.3 
+53.0 + 2.8 


-30.7 + 3.3 
-28.8 + 2.8 


DA 
DA 


9.5 


<0.1% 


1,2 
1,2 



References. — 14Eisenstein et al.l 120061). 2-Kleinman et al., 3 TWegner et al.l l|1987l ). 4 4McCook~ 
IQswalt fc StrunM 1119941) . 6- IGreen et al.l [119861) 

a The secondary in this pair has an unidentified SDSS spectra. 

b J1054+5307 is a wide WDWD previously identified bv lMcCook fc Sionl Jl999l) . 

c The secondary in J1113+3238 is included in the hypervelocity star survey of Lcpinc & Shara 1 120051) . 
^ The secondary in J1412+4216 is included in the hypervelocity star survey of lBrown et al.l <|2007l ). 
c J2115-0741 and J2326-0023 were proposed as wide WDWDs bv lGreavesi (f2005T) . 

spectroscopic coverage in SDSS at only 



Sion (1999), 5- 



confidcnce pairs contains a non-WD. 

As a further test, we calculate the reduced proper mo- 
tion (H g ) for each of our pairs. H g , which combines 
photometric and kinematic information, is an effective 
tool for separating WDs from other objects (|Kilic et al.l 
2006) . Figure [7] shows the reduced proper motions for 
the >99,000 SDSS stars presented in Figure [31 as well as 
lines of constant transverse velocity representing the disk 
and halo WD populations (Vt = 30 and 150 km/s, re- 
spectively). All of the stars in our high-confidence pairs, 
whether spectroscopically confirmed as WDs or lacking 
spectra, lie far from the main-sequence locus and are 
consistent with being WDs. 

We note that two of these pairs, J2115— 0741 and 
J2326— 0023 , were identified as wide WDWDs candidates 
bv lGreavesi (|2005h in a search for CPMBs in SDSS Data 
Release 1. 

In Figure [5J we plot s, the projected orbital separa- 
tion found using the spectroscopically derived distance 
to the primary WDs, as a function of 9 for our high- 
confidence WDWDs. We find no pair with ,s > 10 4 
AU, while our population synthesis predictions are that 
a significant population exists at such separations; we 
also appear only to be identifying very nearby WDWDs. 
However, our meth od cannot yield a complete sample: 
IGirven et al.l (|2011[ ) estimate the completeness of WD 



44% for g < 19 
mag. 

Our tightest pair has 9^8", while we expect binaries 
to exist at smaller 9. This is due to the difficulty in 
matching an object in the SDSS catalog to its USNO- 
B counterpart when the matching radius is of order the 
separation with another object. 

4.2. Comparison with Previously Known WDWDs 

Table [H is a compilation of the properties of all of the 
wide WDWDs reported to date. We include our 12 high- 
confidence WDWDs and 33 systems from the literature 
and present T c s, log g, mass, and cooling age (r), when 
these measurements exist. The available data and their 
quality vary greatly from system to system, but simple 
comparisons can be made between our sample and the 
previously known WDWDs. In Figure El we show /i as 
a function of 9 for the previously known systems and for 
our pairs. Our high-confidence WDWDs are a signifi- 
cant addition to the sample of pairs with small fj, (<200 
mas/yr) and large 9 (>10"). 

Interestingly, four of the systems listed in Tabled] in- 
clude WDs with masses <0.5 Mq. The Galaxy is not 
thought to have had time to produce such low-mass WDs, 
as the youngest WDs in t he oldest Galactic glo bular clus- 
ters have M - 0.5 M Q (|Hansen et al.l [20071. Instead, 
WDs with M < 0.5 M© are likely to form in close bi- 
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TABLE 4 

Wide WDWDs - Derived Quantities 













Prim 


ary 






Secondary 






Name 


s ("; 




ias/yr) a 


T c g (K) 


log g 


M (M ) 


log T 


TaH (K) 


log g M (M Q ) 


log T 


Rof. 









High-C, 


mfrdence W: 


de WDWDs 


Identified in 


This Work b 








J0000— 1051 


16.1 


49.7±8.7 


8598±72 


8.32 + 0.13 


0.80±0.09 


9.23±0.11 










J0029+0015 


8.6 


36.8±6.7 


9947±67 


8.22±0.08 


0.74±0.06 


8.97±0.08 










J0332-0049 c 


18.9 


36.5±9.9 


11012±50 


8.27±0.04 


0.77±0.03 


8.85±0.02 


34288±42 


7.83±0.02 


0.58±0.01 


6.60 + 0.01 


J1002 + 3606 


26.5 


37.2±6.1 


11326±181 


8.05±0.13 


0.63±0.08 


8.68±0.09 










J1054+5307 c 


44.5 


77.6±5.0 


10985±60 


8.08 + 0.05 


0.65±0.03 


8.73±0.04 


11120 


8.01 






J1113+3238 


56.7 


147.3±5.8 


6853±69 


7.56±0.20 


0.38 + 0.08 


8.99±0.08 


7580±88 


8.40±0.15 


0.85±0.10 


9.46±0.09 


J1203+494S 


19.3 


73.6±5.3 


7064±46 


8.05±0.10 


0.62±0.06 


9.21 + 0.08 










J1257+1925 


12.2 


46.9±8.3 


11829±218 


7.72±0.14 


0.46±0.06 


8.46±0.05 


47818±590 


7.85±0.05 


0.62±0.02 


6.33±0.02 


J1412 + 4216 


13.6 


84.4±5.1 


6706±70 


8.76±0.15 


1.07+0.09 


9.64 + 0.01 










J1703+3304 


11.2 


50.9±6.0 


9587±75 


8.11 + 0.11 


0.66±0.07 


8.92±0.09 










J2115-0741 


17.0 


120.7±5.6 


7913±32 


8.16 = 0.05 


0.69±0.03 


9.19±0.04 










J2326-0023 


9.5 


60.2 + 6.1 


7503±67 


8.42±0.12 


0.87±0.07 


9.48±0.05 


10513±46 


8.24±0.05 


0.75±0.04 


8.90±0.04 



sly Identified Wide WDWDs 



LP406-62/63 


28 


480 


5320 


8.0 


0.58 


9.63 


4910 


8.0 


0.58 


9.78 


1,2,4 


LP707-8/9 


12 


172 


















1 


LP647-33/34 


2 


374 


















1 


LP197-5/6 


7 


420 








9.48 








9.94 


1,5 


RE J0317— 853 d 


7 


85 


30000-50000 




>1.1 


~8.4 


16000 


8.19 


0.76-0.84 


~8.4 


6,7 


LP472-70/69 6 


3 




















1 


HS 0507+0434° 


18 


100 


20220 


7.99 


0.62 


7.90 


~12000 


8.1 


0.69 


8.71 


8 


WD 0727+482A/B e 




1340 


5020 


7.92 


0.53 


9.66 


5000 


8.12 


0.66 


9.85 


9 


LP543-33/32 


16 


1800 


4170 


7.65 


0.39 


9.67 


4870 


8.05 


0.6 


9.84 


1,2,4,9 


LP035-288/287 


3 


340 


















1 


PG 0901 + 140 


3.6 


113 


9500 


8.29 




9.10 


8250 








10,11 


PG 0922 + 162 


4.4 


66 


22740 


8.27 


0.79 


7.95 


22130 


8.78 


1.1 


8.41 


10,12 


J0926+1321 d 


4.6 


80 


9500 + 500 




0.62 + 0.10 


8.86 


10482 + 47 


8.54+0.03 


0.79+6 


8.92 


13 


LP462-56A/B 


4 


370 


10240 


8.0 


0.58 


8.90 


8340 


7.5 


0.35 


8.77 


1.2 


LP370-50/51 


13 


215 


















1 


LP549-33/32 


26 


330 


















1,3 


PG 1017+125 


48.8 


30 


















10 


ES0439-162/163 


23 


380 


5810 


8.0 


0.57 


9.52 


4780 


8.0 


0.57 


9.82 


1,2,14 


GD 322 d 


16.1 


87 


14790 


7.87 


0.54 


8.26 


6300 


7.93 


0.54 


9.27 


10,15 


LP322-500A/B 


12 


20 


















1 


LP096-66/65 


18 


300 


















1,16 


L151-81A/B 


2 


50 


14050 


7.96 


0.57 


8.38 


12000 








1,17,18,19 


J1507+5210 d 


5.1 


32 


17622 + 95 


8.13 + 0.02 


0.70 + 0.04 


8.17 


18000 + 1000 




0.99+5 


8.51 


13 


Gr576/577 f 


4.5 


100 


12500 


8.34 


0.8 




9500/8500 




0.39/0.56 




1,20,21,22 


LP567-39/38 


2 


300 


















1 


G206-17/18 


55 


270 


7380 


7.65 




8.93 


6480 


7.75 




9.12 


1.9 


G021-15 f 


58.6 


378 


10000/15000 


8.0/7.4 


0.6/0.35 


8.08 


4750 


8.0 


0.57 




9,10 


GD 392 


45.8 


168 


12220 


9.09 


1.23 


9.22 


~3600 








10,19,23 


G261-43 


1.4 


280 


16000 








5000 








10,24 


HS 2240+1234 


10 


81 


14700 


8.1 






13200 


7.9 






8 


LP701-69/70 


26 


330 








9.34 








9.88 


1,5 


GD 559 


28.7 


134 


















10,25 


LP077-57/56 


3 


409 








9.62 








9.92 


1,5 
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elude a ZZ-Ceti-type variable. 
The primaries in RE .10317-853 and .10926+1321 and secondaries in GD 322 and .11507+5210 ha 
LP472-70/69 lacks a published n, and WD 0727+482A/B lacks a published 0. 
Gr576/577 and G021-15 have both been identified as triple degenerate systems. 



naries whose evolution included a phase of mass trans- 
fer. These four systems are therefore excellent candidate 
triple systems, with the low-mass WD likely to have a 
close-by companion. (Two additional systems in the lit- 
erature are candidate triple systems.) 

The previously known pairs include four WDs that are 
in the Kleinman et al. WD catalog. Of these, only LP 
128-254/255 (J1054+5307) is recovered by our search; 
our primary in this case is LP 128-255. PG 0901+140 
and J1507+5210 were excluded as candidate binaries be- 
cause three of the four WDs in these pairs lack proper 
motion information in SDSS. LP 549—33/32 is not recov- 
ered because the secondary (in our c ase, LP 549—32) has 
colors outside the regi ons defined bylGiryen et al.1 (|201lD 
for DA WDs; indeed, IMcCook fc Sionl (|1999l, classify it 
as a DC WD. The non-detection of three of these four 
systems is therefore not surprising, but confirms that our 
reliance on the Kleinman et al. DA catalog leads us to 
miss a number of wide WDWDs in the SDSS footprint. 



4.2.1. Binary Separation and Stability 

As mentioned in Section [21 wide WDWD orbits are 
thought to be ~5 x larger than those of their progenitor 
systems; mass lost as the stars evolve expands th e or- 
bits from a more compact state (|Greensteinl 1 1 986t ) . We 
compare the projected orbital separations for the WD- 
WDs produced by our population synthesi s to those for 
the w ide main-sequence pairs identified by IDhital et al.l 
(2010) in Figure [TUJ (We also show the individual values 
for our 12 WDWDs; this sample is too small and incom- 
plete for a comparison to the separation distribution to 
be meaningful.) 

The apparent deficit of progenitor pairs for the syn- 
thesized WDWDs with 10 2 < a < 10 3 AU is unsurpris- 
ing, as those progenitor s presumably are tig hter than the 
tightest pairs to which IDhital et al I (|2010f ) are sensitive 
(9 ~ 8"). More interesting is the apparent lack of a signif- 
icant population of WDWDs wider than the widest main- 
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Fig. 8. — Projected separation of our high-confidence WDWDs as 
a function of 9. We use the spectroscopically determined distances 
to the primary WDs to estimate s. Lines of constant distance are 
plotted as dashed lines. 




Fig. 9. — WDWD proper motions as a function of 0. Previously 
known systems are plotted as open stars, while our high-confidence 
pairs are filled stars. The dashed lines indicate the median fi and 
6 values for the previously known WDWDs; our high-confidence 
systems are a significant addition to the sample of pairs with small 
fi and large 6. 



sequence pairs identified by IDhital et al.1 (|2010D . This 
may be because as the widest main-sequence pairs evolve 
to larger separations, they become more likely to be dis- 
rupted by interactions with other bodies in the Galaxy. 

One test of this hypothes is is to compare th e ages of 
our WDWDs to those of the IDhital et al.l (|2010ft binaries. 
We use the cooling ages of our primaries as a rough esti- 
mate of the cooling ages of the systems (we do not know 
whether they evolved first). This cooling age provides 
an upper limit on the WDWDs' lifetimes as truly wide 
binaries, since it is only after both stars have evolved 
into WDs that the pairs reach maximum separations. 
The characteristic lifetimes of the WDWDs in Table 0] 
derived in this ma nner are all >>1 Gyr. By contrast, 
IDhital et all (poll find that many of their pairs have 
characteristic lifetimes <1 Gyr. A larger sample of WD- 
WDs will allow for an improved test of this hypothesis. 

4.2.2. Initial-Final Mass Relation 

While the IFM R is often studied using WDs in 
open clusters (e.g.. [Kalirai et al1l2008f ) or in binaries in 
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Fig. 10. — Normalized orbital separation distributions for the 
WDWDs generated by our population synthesis (das hed line) and 
for th e wide main-sequence binaries identified by Dhi tal et al.l 
1120101 . solid line). For the latter, we use a = 1.1 X s to convert the 
obser ved (projected) or bital separations into true orbital separa- 
tions JDupuy & Liu 2011]). The synthesized systems are those that 
remain after perturbing effects are taken into account (dashed line 
in Figure[T]l. The projected separations for the 12 high-confidence 
WDWDs identified here have also been converted; their orbital 
separations are indicated by the arrows. 

which the WD has a main-sequenc e companion (e.g., 
iCatalan et al.ll2008HZhab et al.ll20l2T ). wide WDWDs can 
also b e used to constrain this relation. iFinlev fe Koesterl 
(|1997f ) estimated the initial mass of the massive sec- 
ondary in the binary PG0922+162 by comparing it to 
similarly massive WDs in open clusters with well-known 
ages. Subtracting the cooling age from the host cluster 
ages gave the pre- WD lifetime of these reference WDs 
and therefore of the secondary in PG 0922+162. This 
was added to the second ary's cooling age t o obta in the 
total age of the system. IFinlev fc Koesterl ([1997) were 
then able to use stellar models to estimate the primary's 
initial mass. The resulting data point in the initial- final 
mass plane has relatively small error bars and can be used 
to discriminate between different theoretical IFMRs. 

However, most of the 45 systems in Table 0] lack spec- 
troscopic information, and many of the WDs with spec- 
tra have large error bars on T e g or log g because these 
were derived from low signal-to-noise spectra. Such er- 
rors can propagate into large uncertainties in the masses 
and cooling ages. While in principle these wide WDWDs 
can be used to constrain the IFMR wit h a method sim- 
ilar to that of IFinlev fc Koesterl (| 19971 ) . we defer a full 
study of this question to a later paper. Spectroscopic 
follow-up of an expanded set of WDWDs, combined with 
ever improving WD atmospheric models, may eventually 
produce a large, homogeneous dataset with which to con- 
strain the IFMR. 

For now, as a simple consistency check, we compare the 
cooling ages and masses of the WDs in the eight systems 
for which both stars have spectroscopically derived T c g- 
and differing masses (and that are not a potential triple 
system) . We expect that in these binaries the more mas- 
sive WD has a larger cooling age. Our simple test holds 
for all eight of these systems. 

5. CONCLUSIONS 
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To expand t he sample of known WDWDs, we adapted 
the method of IDhital et~aH (|2010f) to search the SDSS 
DR7 catalog for WD common proper motion companions 
to spectroscopically confirmed DA WDs out to 9 = 10'. 
We found 41 candidate wide WDWDs. These are pairs in 
which the secondary's colors fall within 0.5<r of the region 
in (g — r) versus (u — g) color space occupied by DA WDs 
(|Girven et al.ll20lTh . We then used two complementary 
techniques to measure the contamination of our sample 
by randomly aligned false matches. We first estimated 
the overall contamination of our sample by false binaries 
by shifting the positions of our primaries several times by 
1° and applying our /i and color criteria to identify false 
companions to the shifted primarie s. We also used the 
Monte Carlo approach developed by IDhital et al.l (|2010p ) 
and searched 10 5 iterations of the line-of-sight to each of 
our primaries for stars whose properties randomly match 
those of these WDs. 

These two tests suggest that, absent other informa- 
tion, the highest probability real pairs in our sample are 
those with 9 < 2' and primaries with random matches 
in fewer than 1% of our rendered lines-of-sight. 13 of 
our candidates meet these criteria; one of these pairs has 
inconsistent spectroscopically derived distance and RV 
measurements for the two WDs, so that our final sam- 
ple includes 12 high-confidence wide WDWDs. (Three of 
these were previously reported as candidate WDWDs.) 

Four of these pairs have SDSS spectra for both DAs, 
wh ile two others have s econdaries classified as DAs 
by iMcCook fc Sionl (|1999f ) . Spectroscopic follow-up is 
clearly needed to confirm the nature of the secondaries 
in the remaining six systems. However, a color selection 
combined with a minimum proper motion requiremen t 
returns a very clean sample of DAs (|Girven "eTaII[20ll ; 
furthermore, based on their reduced proper motions, all 
of our candidates are consistent with being WDs. We 
therefore expect no more than one of our high-confidence 
pairs to contain a non-WD. These systems are a sig- 
nificant addition to the known population with small /i 
( <200 mas/yr) and la rge 9 (>10"). 

iGirven et alj (|2011f ) estimate the completeness of WD 
spectroscopic coverage in SDSS at ~44% for g < 19 mag, 
and it decreases for fainter magnitudes. This is partic- 
ularly unfortunate as the SDSS proper motion catalog 
goes to g ~ 20. Extending the method developed to 



identify the 12 pairs presented here to a photometrically 
selected set of primary WDs could significantly increase 
the number of known WDWDs, which is currently <50. 
An expanded set of WDWDs and spectroscopic follow- 
up may eventually produce a large enough dataset for 
these pairs to realize their full potential as testbeds for 
theories of stellar evolution. 
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